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Abstract 

Using the CLEO II detector at CESR, we have measured the ratio of branch- 
ing fractions B(D+ -r 4>e+v)IB(D+ -r 4>tt+) = 0.54 ± 0.05 ± 0.04. We use 
this measurement to obtain a model dependent estimate of B(Df (^vr"*"). 
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Most measurements of the Df meson branching fractions are normahzed to the clean 
Df — )• (pTT'^ channeL [1] However, the absolute Df — )• ^tt"'" branching fraction is not well 
known, and this limits the precision of these measurements. Here we present a mea- 
surement of Rs = B{Df — )• (pl'^ p) / B (Df — > ^TT"'"), which can be used to extract the 
Df — > (pTT'^ branching fraction by using the known values for the D — )• K*Ip branching 
fractions, the D and Df meson lifetimes and theoretical predictions for the ratio of widths: 
T{Df (Pl+p)/T{D K*l+v). 

The data consist of an integrated luminosity of 1.71 fb~^ of e"'"e~ collisions recorded with 
the CLEO n detector at the Cornell Electron Storage Ring (CESR). A detailed description 
of the CLEO n detector can be found in reference [2]. The data sample contains over two 
million e"'"e~ — )• cc events taken at center-of-mass energies on the T(4S') resonance and in 
the nearby continuum [y/s ~ 10.6 GeV). 

Due to the undetected neutrino, we cannot fully reconstruct Dj" — )• cf)l'^v decays. How- 
ever, there are very few processes which produce both a cf) meson and a lepton contained in 
the same jet. Consequently, this correlation can be used to extract a clean Dj" — )• cf)l'^v sig- 
nal. The backgrounds due to misidentified leptons and from random ^-lepton combinations 
can be reliably estimated, and the possible contamination from other decay modes is shown 
to be negligible. 

We identify cf) candidates by using the decay mode cf) — )• K'^ K~ . In order to suppress 
combinatoric background, cf) candidates are required to have momenta above 1.1 GeV/c. In 
addition, the kaon candidates must have ionization energy loss and time-of-flight consistent 
with that expected for a kaon with the measured momentum. 

The search for leptons is restricted to the kinematic regions in which the lepton identifica- 
tion efficiencies and hadron misidentification rates are well understood. Hence, electron and 
muon candidates are required to be in the fiducial regions | cos 6\ < 0.91 and | cos 6\ < 0.81, 
respectively, where 6 is the polar angle of the track with respect to the beam-axis. In addi- 
tion, electron candidates must have momenta above 0.9 GeV/c and muon candidates above 
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1.4 GeV/c. The only exception is lor muons in the region | cos^| > 0.61 which are required 
to have momenta above 1.9 GeV/c. Electrons are identified by comparing their ionization 
energy loss, time-ol-fiight, and energy deposited in the electromagnetic calorimeter with that 
expected lor an electron with the measured momentum. Electrons Irom photon conversions 
and Dalitz decays ol 7r°'s are rejected by pairing electron candidates with all other oppositely 
charged tracks in the event and rejecting those which have both small separation and parallel 
trajectories at their point of closest approach. Muons are identified by matching charged 
tracks to hits in the muon detectors which lie outside the electromagnetic calorimeter. In 
order to be identified as a muon, a track must penetrate at least 5 interaction lengths of 
steel. Eor leptons in the momentum ranges and fiducial regions considered, the identification 
efficiencies are approximately 92% for electrons and 90% for muons. 

To reduce further the combinatoric background, we require that the (pl'^ momentum 
be greater than 2.4 GeV/c. In order to be consistent with having originated from a Df 
decay, the <f>l'^ candidates must have an invariant mass less than 1.9 GeV/c^. In order to 
suppress the combinatoric background from T(4S') events which tend to be more spherical, 
we require that the ratio of Eox- Wolfram moments [3], i?2 = H2I Hq, is greater than 0.30. 
This eliminates 77% of the T(4S') background whilst retaining 92% of the signal. 

The efficiencies for reconstructing Dj" — )• cf)l'^v decays are obtained from a Monte Carlo 
simulation which takes the predictions of the ISGW model [4] as input. These events are then 
passed through a full simulation of the CLEO II detector and the same event reconstruction 
and analysis chain as the real data. Because of the small value associated with the decay 
cf) — )• K'^ K~ , the kaons tend to overlap in the drift-chamber. This makes it difficult to 
simulate accurately the ionization energy loss measurement. In order to avoid this problem, 
the momentum dependent efficiencies for identifying cf) mesons are obtained from the data 
by comparing the inclusive yield of all ^'s before and after particle identification. These 
efficiencies are then combined with the predicted cf) momentum spectrum from Dj" — )• cf)l'^v 
decays to give the total cf) identification efficiency. EoUowing the above selection criteria, and 
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after correcting for the effects of final-state radiation from the leptons [5], the efficiencies for 
identifying Df — )• (pe'^p decay is fO.5% and for Df — )• (pfi'^p the efficiency is 2.9%. 

Figs f (a) and f (b) show the invariant mass distributions of all K'^ K~ combinations which 
are accompanied by an electron or muon respectively and which pass the above selection cri- 
teria. We fit these distributions with a signal and background function. The signal function 
is a Gaussian function convoluted with a Breit-Wigner function. The background function 
is a phase-space background function [6] which accounts for random K'^ K~ combinations. 
The width of the Breit-Wigner function is fixed to the natural width of the cf) state [8], and 
the mean and sigma of the Gaussian function are fixed to the values extracted from a fit 
to all cf) candidates with momenta above f .f GeV/c. Only the overall normalization of the 
signal function is allowed to vary in the fits. The fits yield 359 it 22 Dj" — )• cf)e'^v and f 23 it f 5 
Dj" — )• (pfi'^p candidates. There are two main sources of background: ^'s accompanied by 
fake leptons [7], and random (pl'^ combinations. 

The background due to fake leptons is estimated by first using the real data to measure 
the momentum dependent probabilities that a hadron will be misidentified as a lepton. These 
probabilities are typically 0.3% for electrons and f .2% for muons. These results are then used 
to randomly label tracks (which do not pass the lepton identification criteria described above) 
as leptons in the data. With this procedure, we extract the number of (pl'^ combinations 
due to misidentified hadrons. For electrons this estimate is 46 it f4 events, while for muons 
it is 27 lb 8. The quoted errors include the contributions from the uncertainties in the 
misidentification probabilities. In order to check these estimates, we examine the invariant 
mass distribution of all (pl'^ candidates. A peak at the Df mass is seen which is due to 
Df — )• (pTT'^ decays in which the pion is misidentified as either an electron or muon. The 
technique of randomly labeled hadronic tracks as leptons yields 7.4 it 2.2 (fnr events where 
the TT track is misidentified as a lepton. This is in good agreement with the 4.5 it 3.6 events 
found in the lepton sample. 

For the range of lepton momenta considered, random (pl~^ combinations come from two 
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sources: from e~'~e~ — )• cc events in which a ^ is produced in the fragmentation process and is 
combined with a lepton from the semileptonic decay of the charmed hadron in the same jet, 
and from T(4S') decays in which a ^ is produced in the decay chain of one of the B mesons 
and is combined with a lepton from the semileptonic decay of the other B meson. A (f> and 
lepton which originate from the decay of the same B meson do not contribute since they 
tend to be emitted back-to-back and thus have too large an invariant mass. The background 
from random (pl'^ combinations is estimated using the Monte Carlo simulation. However, 
this is complicated by the fact that the (f> production rate from both fragmentation and B 
meson decays is not well known. For this reason, an attempt is made to scale the Monte 
Carlo prediction to account for the (f> production rate observed in the data. 

In the continuum, not only the rate of (f> production, but also the correlation between the (f> 
and the charmed hadron direction is important. The agreement between the data and Monte 
Carlo is investigated by considering how often a ^ is produced in the same hemisphere as a 
fully reconstructed D meson. Both _D° and D*"*" mesons are considered, and are reconstructed 
using the following decay chains: _D° — )• K~Tr'^ and D*'^ — )• _D°7r"'"; _D° — )• K~Tr'^. The 
reconstructed D mesons are required to have momenta above 2.5 GeV/c in order to account 
approximately for the range of D momenta which are expected to contribute leptons in the 
momentum range of interest. For this particular study, the (f> momentum criterion is relaxed 
to 0.8 GeV/c in order to provide sufficient statistics. In doing this we have assumed that 
the (f> momentum distribution is well reproduced by the Monte Carlo, and that it is the rate 
of (f> production in the fragmentation process which contributes the greatest uncertainty. 
Both the number of D mesons and the number ^'s are obtained by fitting their invariant 
mass distributions. False combinations due to the D meson combinatoric backgrounds are 
accounted for by subtracting the number of ^'s found when using the D mesons invariant 
mass sidebands. In the real data 0.17 it 0.11 ^'s are found for every 1000 reconstructed D 
mesons. This is to be compared with 0.16 it 0.02 for the e"'"e~ — )• cc Monte Carlo. The ratio 
of these two numbers is 1.0 it 0.7, so no correction is applied in this case. The simulation 



6 



predicts a background of 12 it 8 and 1.8 it 1.2 events for electrons and muons respectively, 
where the errors include the uncertainty in the above ratio. 

The background from random (pl'^ combinations in T(4S') decays is estimated in a similar 
manner. In this case the directions of the (f> and lepton are uncorrelated. For this reason it is 
sufficient to compare the number of ^'s with momentum above 1.1 GeV/c in the continuum 
subtracted T(4S') data with that observed in the T(4S') BB Monte Carlo. In the real data 
5.0 lb 0.5 ^'s are found per 1000 BB events to be compared with 5.3 it 0.1 in the Monte Carlo. 
This gives a correction factor of 0.95 it 0.08. After applying this correction, the predicted 
background is 19 it 2 events for electrons and 9 it 1 events for muons. 

Figs 2(a) and 2(b) show the number of ^'s which fall in each (pl'^ invariant mass bin for 
electrons and muons respectively; where the number of ^'s has been extracted from fits to the 
K'^ K~ invariant mass distributions. The combined background estimates are also shown, 
as well as the simulated predictions for the signal shapes which have been normalized to the 
number of candidates extracted from the fits to the K~ invariant mass spectra. It can be 
seen that the predicted signal shapes are in good agreement with the data. The background 
estimates can be checked by comparing the predicted number of candidates which fall outside 
of the signal region with the number actually observed. For electrons we predict 8 it 1 events 
in the region 2.0 < Af^;+ < 3.5 GeV/c^ and observe 12 it 7, and for muons we predict 7 it 1 
and observe 8 it 5. Both predictions are in good agreement with the data. 

We have also estimated the possible contamination from the decays D'^ — )• cf)K^l'^v, 
Df — )• (prjl'^p and Df — )• (pTnrl'^p. The decay Df — )• cf)TT^l'^v is forbidden from conservation 
of isospin. The Feynman diagrams for the first two processes are shown in Figs. 3(a) and 
3(b), respectively. For the first decay, we first estimate an upper limit on the number of 
D'^ — )• ^A^°/"'"z/ decays in the data by making the conservative assumption: B{D'^ — )• cf)K^l'^v) 
IB{D^ K^l^v) ~ a ■ B{D° {K*tt)~ ji^v) IB{D^ K~ ji^v), where a is a suppression 
factor because in the first numerator it is an ss pair which must be popped from the vacuum 
as opposed to a light-quark pair [9]. This assumption is motivated by the similarity of 
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the decay diagrams. The latter ratio has been measured by the E653 collaboration [10], 
B{D° {K*7r)-fi+iy)/B{D° R-p+v) < 0.04 at the 90% confidence level. Taking the 
number of D'^ — )• K^l'^v events in our data sample to be 60,000 [11], assuming a = 1/3, 
and including the simulated acceptance for D'^ — )• cf)K^l'^v decays, we estimate less than 
one background event from this source. The contribution from Dj" — )• (prjl'^p decays is 
estimated in a similar manner. Here we make use of the same E653 result and assume: 
BiDf (t>f]l+v)IB{D+ (t>l+v) c^a- 13- B{D'' (A^*7r)-/i+//)/B(i)° ^ K-iJ+v), where /3 
accounts for OZI suppression [12] of the first numerator. Making the assumption fi = 1/10, 
and including the simulated acceptance for Dj" — )• (prjl'^p decays, we again estimate much less 
than one background event from this source. The contribution from Df — )• (pTnrl'^p decays 
should also be very small for similar reasons. Therefore, it is assumed that the background 
from these decay modes is negligible. 

After subtracting all backgrounds, we find 282 it 22 Df — )• (pe'^p and 85 it 15 Df — )• (pfi'^p 
candidates which fall in the Df signal region M^i+ < 1.9 GeV/c^. After correcting for 
the detection efficiencies in each channel and for the (f> — )• K'^ K~ branching fraction [8], 
the efficiency corrected yields are 5460 it 430 for electrons and 6000 it 1000 for muons. A 
breakdown of the yields in each channel is given in Table I. To combine these two numbers, 
we take a weighted average, after first correcting for the fact that the muon rate is predicted 
to be 5% lower than that for electrons because of the reduced phase-space [13]. Therefore, 
our result is given in terms of the effective yield in the electron channel which is 5580 it 400 
events. 

In order to limit the systematic effects which stem from the selection criteria, the number 
of Dj" — )• cf)TT'^ decays is measured in a similar manner. Again, the cf) candidates are required 
to have momenta above 1.1 GeV/c . To account approximately for the fact that no neutrino 
is produced in this decay, we require the ^tt"'" momentum to be greater than 2.7 GeV/c [14]. 
We then require the K'^ K~tt'^ invariant mass to be within ±25 MeV/c^ of the known Dj" 
mass [8]. The efficiency for detecting — )• ^tt"*" decays following these selection criteria is 
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17.4%. 



The result of the fit to the K'^ K~ invariant mass distribution is shown in Fig. 4. We find 
1049 lb 35 candidates. Also shown is the result of the fit to K'^ K~ combinations from the Dj" 
mass sidebands, which is used to estimate the contribution from random ^tt"'" combinations. 
We find 163 it 18 candidates due to these random combinations. After subtracting this 
background, the efficiency corrected yield is 10, 370 it 460 events. 

Finally, since we have already corrected for efficiencies, the ratio of branching fractions 

is. 



where the first error is statistical, and the second is an estimate of possible systematic 
effects. This systematic error includes: the uncertainty in the number of fake leptons (6.3%), 
the uncertainty in the level of continuum charm background (2.7%), the uncertainty in 
the level of BB background (0.9%), the uncertainty in the lepton identification efficiency 
(2.5%), the uncertainty in the cf) identification efficiency (1.0%) and that due to the limited 
number of Monte Carlo events which were used for the efficiency estimates (2.7%). We 
have also considered our sensitivity to the Dj" production mechanism by using the predicted 
efficiencies for Dj" mesons produced in D*'^ decays. However, since both the Dj" — )• cf)l'^v 
and the Dj" — )• ^tt"'" efficiencies are affected in the same manner, the effect is small (0.5%). 
For this analysis the ISGW model was used to generate semileptonic decays in the Monte 
Carlo simulation. The uncertainty associated with this choice of model was investigated 
by adopting the iV — A) prediction for the lepton momentum spectra. From the resulting 
change in the predicted Dj" — )• cf)l'^v acceptance, we assign a systematic error of 1.8%. 
Various background shapes have also been used to extract the number of cf) mesons. In all 
cases Rg changes by less than 1.1%, which is taken to be the systematic error. After adding 
these estimates in quadrature, the total systematic error is 7.6%. Table II compares this 
result with those of previous measurements [15-17]. 

Having measured i?^, we extract the Dj" — )• ^tt"'" branching fraction by using the theo- 




10,370 ±460 



5580 ± 400 



0.54 ±0.05 ±0.04, 
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retical value lor Ts = ^'{Df (/)l~^u)/T{D K*l^v). We can write 



V) X 



where Tj^+ and T£,o are the Dj" and _D° lifetimes, respectively. JF^ is predicted by various 
quark models [18,19]. Here we choose to adopt the result ol the modified ISGW model 
J^s = 1-00, since to date this is the only model which can account for the measured value for 
B{D K*l+v)IB{D Kl+v) [19,20]. Past experiments have used J^, = 0.9 for this ratio. 
We use CLEO II measurements for all quantities except the _D° and Dj" lifetimes. In this way 
some of the systematic errors cancel and the problems associated with averaging the results 
of many different experiments are avoided. We found B{D^ — )• K*~ e'^ v) j B{D^ — )• K~tt'^) = 

0. 61 lb 0.07 [11], and when this is combined with our measurement, B{D^ — )• K~tt'^) = 
(3.91 ±0.19)% [21], and with the E687 measurements, r^o = (4.13 ±0.05) x lO'^^ g [22] and 

= (4.75 ± 0.21) X 10-^3 g [23], we obtain B{D+ c^e+v) = (2.74 ± 0.36)%. Using our 
measurement of i?^, we obtain: B{Df — )• ^tt"'") = (5.1 ± 0.4 ± 0.4 ± 0.7)%, where the first 
error is from the statistical error on i?^, the second from the systematic error on i?^, and 
the third from the uncertainty in the Df — )• (pe'^p branching fraction. This measurement 
is greater than but consistent with previous estimates and upper limits on B{Df — )• ^tt"'") 
[8,24]. 

In conclusion, we have measured B{Df — )• (pe'^ p) / B (Df — )• ^tt"'") = 0.54 ±0.05 ±0.04. By 
using the theoretical prediction J^, = 1.00, we find B{Dj (f)Tr+)= (5.1 ± 0.4 ± 0.4 ± 0.7)%. 
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FIGURES 
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FIG. 1. Fits to the K'^K~ invariant mass distributions for (a) K'^K~e'^ and (b) K'^ K~ ji 
combinations which lie in the _D+ signal region Mj^+j^-i+ < 1.9 GeV/c^. 



13 



50 



1 1 1 1 1 1 1 1 1 1 1 1 r 



(b) 




0,990 1.010 1.050 1.050 

K"K" Moss (GeV/c^) 



14 



100 




1.0 



1.5 2.0 2.5 

r Mass (GeV/c^) 



3.0 



3.5 



FIG. 2. Invariant mass of (a) (^e"*" and (b) (^^u"*" combinations. The data points are obtained 
by fitting the K'^K~ invariant mass distributions for each cj)!'^ invariant mass bin. The solid 
histograms show the sums of the predicted backgrounds and the simulated signal shapes. The 
dashed histograms show the background contributions. 
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FIG. 3. Feynman diagrams for the possible background modes, (a) (j)K^l'^v, in which 

an ss pair must be "popped" from the vacuum, and (b) _D+ (j)ril'^v, which is OZI suppressed. 
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FIG. 4. Fit to the K'^K~ invariant mass distribution for K'^K~'it'^ combinations which lie 
within ±25 MeV/c^ of the Df mass. The dashed histogram shows the contribution from the Df 
mass sidebands. 
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TABLES 



TABLE L Summary of yields. The errors quoted in this table are statistical only. 



Decay mode 


D+ -r 4>e+v 




Total candidates 


359 ± 22 


123 ± 15 


Lake lepton background 


46 ± 0.3 


27 ±0.8 


Continuum cc background 


12 ± 0.4 


1.8 ± 0.1 


BB background 


19 ± 0.8 


9± 0.5 


Background subtracted 


282 ± 22 


85 ± 15 


Efficiency, e ■ B {%) 


5.16 


1.42 


Efficiency corrected yield 


5460 ± 430 


6000 ± 1000 



TABLE IL Comparison of this result with those of previous experiments. We have increased 
the E687 result by 5% since only muons were used in their analysis. 



Experiment 


Events 


Rs 


CLEO 1.5 [15] 


54 


o.49± o.iotEJ:^!^ 


ARGUS [16] 


104 


0.57 ±0.15 ±0.15 


E687 [17] 


97 


0.61 ±0.18 ±0.07 


This result 


367 


0.54 ±0.05 ±0.04 
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